Optimal quantum circuits for general phase estimation 
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We address the problem of estimating the phase (j> given copies of the phase rotation gate u^. 
We consider, for the first time, the optimization of the general case where the circuit consists of 
an arbitrary input state, followed by any arrangement of the A'^ phase rotations interspersed with 
arbitrary quantum operations, and ending with a POVM. Using the polynomial method, we show 
that, in all cases where the measure of quality of the estimate 4> for (j> depends only on the difference 
(j) — 4>, the optimal scheme has a very simple fixed form. This implies that an optimal general phase 
estimation procedure can be found by just optimizing the amplitudes of the initial state. 

PACS numbers; 03.67.-a 



The possibility of encoding information into the rela- 
tive phase of quantum systems is often exploited in sev- 
eral quantum information processing tasks and several 
kinds of applications. For example, it was shown that in 
most existing quantum algorithms the information to be 
retrieved after the computation is contained in relative 
phases [Q. Moreover, information is encoded into phase 
properties in some quantum cryptographic protocols [^, 
and in some precision measurements, such as the schemes 
on which atomic clocks are based Therefore, the is- 
sue of estimating the phase in the most efficient way is 
of great interest. 

We phrase the phase estimation problem as follows. 
Let Urf, be a single qubit gate that, in a prescribed 'com- 
putational' basis {|0), |1)} maps the state |0) to |0) and 
|1) to e"^|l). We assume we have no prior knowledge 
about <j). The objective is to estimate (j) using some pro- 
cedure that will output some guess (p. We characterize 
the quality of a estimate by a "cost function" C((/), 0), 
which specifies the penalty associated with guessing (p 
when the actual phase is p. We are given N identical 
single qubit quantum gates u^, and the goal is to use 
these gates along with any other operations in order to 
produce an estimate of (p. The optimal procedure is the 
one that has the minimum expected cost. 

Most of the previous work on phase estimation assumes 
some fixed state encoding the phases, and the only thing 
to be optimized is the final POVM measurement ^, ||, 
[iot . More recent work fixes the way the phase gates 
are applied and optimizes the choice of input state and 
final POVM. 

The crucial point is that in this paper we are not re- 
stricted to preparing some input state, then applying all 
of the phase rotations, and then performing an optimal 
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POVM. We consider, for the first time, the case where 
one has full freedom over how to use the phase rotation 
gates in an experiment designed to optimally estimate 
the phase. Any realistic experiment of this type can be 
viewed as computation, completely specified by a quan- 
tum circuit acting on some finite number of qubits and in- 
volving, apart from the N copies of the gates, some fi- 
nite number of arbitrary quantum gates of our choice, fn 
fact, many quantum algorithms, including Shor's quan- 
tum algorithm for factoring integers, can be phrased in 
terms of such phase estimations This has origi- 

nally provided motivation for this work. 

We assume the phase p is chosen uniformly from [0, 27r) 
|]l8| and that a suitable quantum circuit containing N 
copies of the gates outputs some value y with prob- 
ability Pr(j/|0). From y we infer, following a prescribed 
rule, the estimate (py. The quality of the whole procedure 
is quantified by the expected cost C, given by 

The next part of this paper describes a very simple 
procedure for estimating p, that only requires one to op- 
timize the choice of initial state to an otherwise fixed 
procedure. The rest of the paper then reduces the very 
general case we have described above to this very simple 
case. 

For a given cost function, the quality of an estimation 
procedure depends on both Vv(ii\<p) and the inference rule 
y ^ Py The optimal protocol gives the minimum possi- 
ble average C. We restrict attention to cost functions C 
that depend only on (p — py, and therefore adopt the no- 
tation C{p, (py) ~ C{p — (f)y). We will make only the fol- 
lowing very weak assumption on the cost function (which 
corresponds to a more general class of cost functions than 
the "Holevo" class that is typically considered j^) 

/ d(t)\C{P)\ < oo. (2) 
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FIG. 1: A simple approach for estimating (p, in the case that 
N = 7. Optimizing over the input amplitudes aj, produces 
an optimal estimate of 0. 



FIG. 2: Most general approach for estimating the unknown 
phase rotation. This subsumes Procedure |l|, as well as more 
complicated procedures with classical feedback. 



We will deal with specific cost functions later on. 

Let us start by describing a simple and natural ap- 
proach for estimating the phase 0, illustrated in Figure 

Procedure 1 • Prepare m qubits in state \x) = 

J2j=o'^j\3)' ^ — 2™ — 1. The exact values 
of aj depend on the cost function to be maximized. 

• Apply the gates to effect 

N N 

i=o j=o 

• Apply the inverse quantum Fourier transform to ob- 
tain 

2"-i I N 2 \ 

2-W2^ [Y^a,e^^^^-^A\y)^ (4) 

y=0 \j=0 j 

measure y and calculate the estimate (i>y = 

The surprising claim is the following. Given any func- 
tion C satisfying Eq. (Q), the minimum of C obtained 
by optimizing the aj in Procedure |l| is the infimum of 
all values obtainable by any realistic experiment (as we 
described above and illustrate in Figure ^ . It is impor- 
tant to also note that apart from the preparation of the 
initial state, the above procedure can be implemented 
using the N black boxes and a number of elemen- 
tary gates polynomial in log(A^). Efficient preparation 
of states is discussed in |l2|. Exact implementation of 
quantum Fourier transforms is discussed in [ p^, a. nd ar- 
bitrarily good approximations are discussed in |l^ . 

The remainder of this paper will prove this claim by a 
sequence of reductions. 

Let us start with a very general circuit (Figure H) which 
uses m+d qubits, where m and d can be arbitrarily large. 
The first m qubits are measured after the computation, 
yielding the output < y < 2'" — 1, whereas the remain- 
ing d qubits are discarded. 

Since we are allowing arbitrarily many extra "ancilla" 
qubits, and since any classical feedback scheme can in 
principle be implemented by a unitary operation using 



a sufficiently large ancilla, the family of schemes that 
can be implemented by a quantum circuit of this form 
includes any scheme using finite dimensional state spaces. 

For convenience, we let the output y correspond to the 
phase estimate = for some M, M < 2™. For 
a fixed approximation scheme (using finite means) and 
cost function satisfying Eq. (|^) , and assuming a uniform 
prior distribution of the 0, this simplifying assumption 
will give us a scheme with expected cost that is at most 
C + eM, where em ^ as Af ^ oo, and C is the lowest 
expected cost for any possible scheme. Thus the infi- 
mum of the C over all such restricted schemes equals the 
infimum of the C over all possible such schemes. 

In fact we will also show later that as long as M > 

-|- 1 , this simplifying assumption does not cost us any- 
thing. That is, the infimum of the expected costs of all 
the schemes using the inference rule y t-^ "^jj- for any 
M > TV -I- 1 is the infimum of the possible expected costs 
using any inference rule y ^-^ (j)y. 

Suppose we came up with a general circuit that per- 
forms an estimation of (j) according to some prescribed 
set of criteria. Let us first show that such a circuit is 
equivalent, for our purposes, to another one, which has 
much simpler structure. 

The state at the output of the circuit can be written 

as 

M 2'^-l 

a{y,z,(j))\y)\z), (5) 

j,=0 2 = 

where each amplitude a.{y^ z, cj)) is a polynomial in e"^ of 
degree at most N 

a{y,zA) = j:^^^e-t (6) 

This fact follows just as in [Q where the polynomial 
method is applied to an oracle revealing one of many 
Boolean variables. This paper is the first time the poly- 
nomial method is applied to oracles with continuous vari- 
ables. Since we assume that the cost function is of the 
form C{(j),(j)y) = C{(j) ~ 4'y)j then we can also assume, 
without loss of generality that the optimal conditional 



probability 



Pr(2/|(/.)=^|a(y,z,0)P (7) 

z 

and therefore 



depends only on the difference <// — -j^ 
equals 



Pr(0|(/. - 27ry/M) = ^ |«(0, z> - 27ry/M)\\ (8) 



To simplify the notation, we let aj{z) — aj{0, z). There- 
fore, a circuit that produces amplitudes 



1 ^ 

\a{y,z,(l)) \ = j^\'^a:3iz)e' 



(9) 



3=0 




also leads to the optimal estimation of (f>. Thus the fol- 
lowing simple estimation procedure, whose circuit is il- 
lustrated in Figure |^, performs equally well: 

• Prepare m + d qubits in state \x) = 

Z]z=o^ ZiLo "i(^)li)l^)- ^'^^ ^^^^ preparation 

to be possible m has to be chosen such that 
N < 2™. 

• Apply the gates to effect on the first m qubits 

2''-! N 2"^-! N 

z=0 j=0 z=0 ■j=0 



Apply the inverse quantum Fourier transform, ^ 
\j) ^ (l/VM)Eyio'e-'^|2;), to obtain 

E E E",We'^'^^-^^ ) |y>k>, (11) 

z=0 y^O \i=0 



and measure y. 

The following two observations lead to further simplifi- 
cations. 

Let us first notice that in this procedure the role of 
the d auxiliary qubits is restricted to the initial prepara- 
tion of the most general state of the first m qubits (all 
subsequent operations are restricted to these m qubits). 
Usually, such a state is described by a density operator 
which can always be expressed as a mixture of pure states 
of the first m qubits. The average cost in this case is given 
by the average, over the mixture, of individual expected 
costs pertaining to the pure states in the mixture. Thus 
the expected cost for the mixture cannot be greater than 
all of the individual expected costs for the contributing 
pure states; hence either some of the contributing costs 
are greater or they are all equal. In either case a judicial 
choice of a pure state of the m qubits leads to equally 
good or better phase estimation. This argument implies 
that without loss of generality we can restrict our circuit 



FIG. 3: Without loss of generality, we can assume our phase 
estimation procedure has the form illustrated in this figure, 
where 2™ > M > A'' + 1. The top register contains m qubits, 
and the bottom d qubits are ancilla qubits that may be en- 
tangled with the first m qubits, but are discarded. 



to only m qubits (plus some ancilla bits used to imple- 
ment using N copies of u^) and run the estimation 
on pure states. 

Secondly, in the description above the quantum Fourier 
transform is parameterized by A/, where -t- 1 < M < 
2", but in fact any Af > iV + 1, in particular M = 2™, 
will work equally well. To see this consider a cost function 
of the form C(0 — (j)y). The expected cost, for some Af, 
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where 4>' = 4> — . The expected cost does not depend 
on M as long as M >N + \. 

Recall that we mentioned in the introduction that as 
Af — > oo the difference between the optimal C assuming 
the inference rule y ^ (jiy = and the optimal C 
without such an assumption is e^i 0. Since we have 
just shown that for all Af > -I- 1, the expected cost C 
is constant, this means that the difference eM is in fact 
once Af > A'' -I- 1. In other words, assuming 4>y = "^jj- 
does not cost us anything as long as we use M > N + \. 

It is clear that the exact value of the expected cost 
now depends only on the aj values, that is, on the initial 
state, which means that given a specific cost function all 
we have to do is to choose an optimal initial state. 

We emphasise that the schemes in Figs. ^ and || pro- 
vide an optimal covariant estimation scheme even for gen- 
eral cost functions not necessarily of the Holevo class. 
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Indeed, in ||l7|] it is proved that there exists a discrete 
POVM achieving the same average cost of any covari- 
ant continuous POVM. By the Naimark theorem Q this 
means that there exists an orthogonal measurement on 
the system and an ancilla achieving the discrete POVM, 
and a further extension aUows to have the POVM as 
rank-one (projectors of rank r > 1 can be easily con- 
nected to rank-one projectors via a suitable controUed- 
unitary interaction with an ancilla). 

Let us now address the problem of optimal input 
states for two different cost functions. First we look 
at the minimization of the "1— Fidelity" cost function 
CF{(j),(f)y) = sin^[((/) — 0y)/2]. The minimum cost is 
achieved with the initial state 

N 

i-T""^') = E\/S--(^)b-)- (12) 

The error in fidelity of this protocol goes to zero ac- 
cording to the square of the number of black boxes used 
Cf = 0(1/N'^). It is interesting to note that the fidelity 
of the more conventional approach to phase-rotation es- 
timation with the uniform initial state (a,- — ,}, , ^ 
for all j) only tends to zero linearly in N. That is 
CF=n(l/N). 

Another cost function that is commonly used is the 
window function that allows any error smaller than 5: 
^) = if < 5, but Cfyicj), 4>) = 1 if 10- ^1 > 

5. The minimisation of this cost leads to optimal states 



with amplitudes aj = + 1, which corresponds to 

what is effectively used by Shor's algorithm j!], |l^ 
and provides an expected cost in 0{j^). 

In this paper we have addressed the general problem 
of finding the optimal estimating procedure for the real 
parameter given N copies of the single qubit phase 
rotation within a general quantum circuit in finite 
dimensions. We considered the general case where the 
circuit consists of an arbitrary input state followed by 
any arrangement of the N phase rotations interspersed 
with arbitrary quantum operations. The main result was 
the proof that in all cases, and for any covariant cost 
function, we want to use the optimal phase estimation 
procedure is equivalent to a quantum Fourier transform 
in an appropriate basis. 

Our result is very general, and gives a recipe for finding 
the best achievable phase estimation for a given cost func- 
tion. In practice, once we know the minimum cost possi- 
ble, one can also search for and use easier-to-implement 
phase estimation procedures that achieve the same, or 
similar expected cost. Due to the generality of our main 
result, it will surely find many other interesting applica- 
tions in physical and computational scenarios. 

This is the first application of the polynomial method 
to "black-boxes" encoding continuous variables, in this 
case, one real parameter. The method can also be applied 
to several real parameters, as well as combinations of real 
and discrete parameters. 
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